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Abstract 

The limits of the Second Law of Thermodynamics, which reigns undisputed in the macroscopic 
world, are investigated at the mesoscopic level, corresponding to spatial dimensions of a few mi- 
crons. An extremely simple isolated system, modeled after Callen's adiabatic piston, [J can, under 
appropriate conditions, be described as a self-organizing Brownian motor and shown to exhibit a 
perpetuum mobile behavior. 

PACS numbers: 05.70.Ln, 05.40.Jc 
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I. INTRODUCTION 



As first clearly stated by Schroedinger "a living organism tends to approach the danger- 
ous state of maximum entropy, which is death. It can only keep aloof from it, i.e., alive, by 
continually drawing from its environment negative entropy. What an organism fed upon is 
negative entropy" . [2| Is it possible for an organism immersed in a thermal bath to be insu- 
lated, and thus avoid thermalization, and still be able to reduce its entropy? Unfortunately, 
the above scenario contradicts one of the most cherished laws of physics, that is the Second 
Law of Thermodynamics, which has victoriously resisted all the attempts aimed at finding 
a particular case where it could be violated. However, one has to consider that the Second 
Law has been formulated in the context of macroscopic physics and it is in this context 
that it has been successfully applied and verified. While the Second Law, being inherently 
statistical in its nature, cannot carry over to microscopic cases where the number of involved 
particles is too small, its limits of validity are not well understood in the transition region 
bridging the macroscopic to the microscopic world, where the system dimensions are small 
but still the number of particles large enough to justify the use of Thermodynamics. 

In particular, some intermediate mesoscopic regime, like, e.g., the biological one associ- 
ated with typical cell dimensions (about one micron), is still terra incognita and open to 
possible surprises. Recent years have indeed witnessed a growing interest in the Thermody- 
namics of small-scale non-equilibrium devices, especially in connection with the operation 
and the efficiency of Brownian motors. t 3| These are devices aimed at extracting useful work 
out of thermal noise and are the descendents of the famous ratchet mechanism popularized 

n 

by Feynman |4J in order to show the impossibility of violating the Second Law of Ther- 
modynamics. From a quantitative point of view, they are usually described in terms of a 
variable x(t), which can be identified as the trajectory of a "Brownian particle" of given 
mass M obeying Newton's equation of motion. The forces acting on the "particle" are those 
resulting from some kind of prescribed potential V(x), from the viscous-friction term — , yx(t) 
and from the randomly fluctuating force (Langevin's force) associated with thermal noise. 
Even in the presence of a spatially asymmetric (ratchet-type) potential V(x), no preferential 
direction of motion is possible and no perpetuum mobile mechanism leading to a violation of 
the second law. 4] In this paper, we show that a very fundamental thermodynamic system, 
the so-called "adiabatic piston", first introduced by Callen,jl| can be actually described as 
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FIG. 1: The adiabatic piston: an insulating cylinder divided into two regions by a movable, 
frictionless and insulating piston. 

a self-organizing Brownian motor. Under specific conditions, the system can be analytically 
modeled as a harmonic oscillator undergoing Brownian motion. The relevant variable of the 
system is the position x(t) of the moving wall with respect to its central position L/2, the 
harmonic potential being not prescribed but naturally emerging, as well as the friction term, 
as the result of the internal dynamics of the system. 

We recall that the adiabatic-piston problem deals with a system consisting (see Fig.l) of 
an adiabatic cylinder of length L, divided into two sections A and B, containing the same 
number n of gas moles in different equilibrium states, by a fixed frictionless adiabatic piston 
. At a given instant t=0, the piston is let free to move and the problem is to predict the 
final equilibrium state of the system. This deceivingly simple question, which cannot be 
actually answered in the frame of equilibrium thermodynamics, has been firstly solved in 
1996 in the frame of a simple kinetic model and has since then given rise to a number 
of papers of increasing sophistication. Both analytic [fj, Q, 0, E| and molecular dynamics 



approaches 
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131 ] reveal the existence of two successive stages of evolution, the 
first one leading the system to a state of mechanical equilibrium (where the pressures in 
both sections are equal) and the second, intrinsically random and occurring over a much 
longer time scale, to a final state of both mechanical and thermal equilibrium (where the 
average position of the piston is equal to L/2 and the average gas temperatures in the two 
sections become equal). The above conclusions have been basically drawn by investigating 
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the average behavior of the significant variables describing the system evolution but less 
attention has been devoted to the fluctuations around these average quantities, which are 
somehow assumed to be small. By extending and clarifying the results of a previous paper 
Q|, we wish to show how this is not always the case. More precisely, if one starts from 
specific initial conditions, then for spatial dimensions of the cylinder corresponding to the 
mesoscopic range (about 1 micron) and for values of the ratio \x = M/M g between the piston 
and the gas mass lesser than one, fluctuations around equilibrium position can be a sizeable 



fraction of the average quantities, larger by a factor y Mj m (m being the mass of the single 
gas molecule) than the ones expected if the system were in thermal equilibrium. This makes 
rather questionable the statement that in this case a final equilibrium configuration can be 
actually reached. Actually, the evolution of the system from an initial state in which both 
x(t) = X(t) — L/2 and x(t) are zero and the temperatures are equal in the two sections 
exhibits unexpected features: the system fails to settle down to a well-defined equilibrium 
state since the piston never comes to a stop but keeps wandering symmetrically around the 
initial position, performing oscillations whose mean-square value (evaluated over an ensem- 
ble of many replica) is much larger than that pertaining to standard thermal fluctuations. 
According to this behavior, the system appears to deserve the name of perpetuum mobile, 
even if there is no preferential direction of motion ((x) = (x) = 0). This unusual behavior 
can be associated with a possible challenge to the Second Law. In fact, a suitable process 
can be conceived through which the entropy of the isolated adiabatic piston system turns 
out to undergo a significant decrease (see Sections IV- V). In practice, as we shall see below, 
no significant entropy decrease occurs outside a mesoscopic regime corresponding to spatial 
dimensions of a few microns; in fact, the time interval t as over which the significant entropy 
decrease occurs depends on the fourth power of L (t as cx L 4 ): thus, any deviation from the 
above spatial scale gives rise to such large t as as to affect the validity of our model. We note 
that the coincidence between the spatial scales over which the model is valid and a typical 
biological dimension is definitely intriguing. 

II. THE STOCHASTIC ADIABATIC PISTON 

The investigations of the adiabatic-piston dynamics appeared in the last ten years (see, 
e -g-> P>y> Q fl El) have placed into evidence two distinct stages of temporal evolution. The 
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first one occurs, on a relatively short time scale, when the system is let free to evolve from 
an initial configuration in which the two sections are in thermal equilibrium at different 
pressures, and leads to a final mechanical equilibrium (equal pressures) configuration that 
depends on the initial conditions. The second one is characterized by a stochastic dynamic 
evolution during which the common pressure on the two sides remains essentially constant, 
and may eventually lead to a final equilibrium configuration in which also the temperatures 
attain a common value and the piston is in the central position. In this section, we confine 
ourselves to the second stage by investigating the piston dynamics starting from a specific 
initial configuration characterized by x(t = 0) = x(t = 0) = and by a common temperature 
T in the two sections. In particular, we wish to focus our attention on the fluctuations of 
the piston positions around L/2. In effect, while (x(t)) and (x(t)) turn out to vanish at all 
times, as a priori dictated by symmetry, there is a range of values of the system parameters, 
corresponding to a mesoscopic regime where the mean-square value {x 2 (t)) l l 2 can be a not 
negligible fraction of the piston half-length L/2, a quite surprising result which seems to 
indicate that no effective final equilibrium position can be reached. 

We start from the equation describing the deterministic piston motion in the presence of 
a finite pressure difference, that is [5( 



nR (T A T B \ l8nRM g / ^T A ^T b \ M g ( 1 1 x • 2 



where R,M g and M are the molar gas constant, the common values of the gas masses in A 
and B, and the piston mass, respectively, while the dot stands for derivation with respect 
to time. It has been derived by assuming that the velocity distributions of the ideal gas 
in the two sections are Maxwellian, i.e., that thermal equilibrium is continuously restored 
on a time scale much shorter than the characteristic time scale of the piston motion, an 
assumption which has to be checked a posteriori. We observe that Eq.(l) is identical to 
the one obtained by Munakata and Ogawa (Eq. (22)of Ref. Q])but for the extra-term 
containing X 2 . This term, which is negligible for finite pressure-difference in the two sides, 
turns out to play a fundamental role when studying the piston dynamics for equal pressures, 
i.e., in the situation investigated in the present paper. In this case, Eq.((TJ) has to be suitably 
modified. In effect, besides setting equal to zero the first term on its R.H.S. associated 
with pressure difference, we have to explicitly take into account the random nature of the 
hits suffered by the piston because of the collisions with the gas molecules. To this aim, 



we resort to Langevin's approach (see, e.g. based on the introduction of an ad hoc 

Langevin acceleration a(t), so that Eq.Q is superseded by 
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where is the Boltzmann constant and N = M g /m is the number of molecules of mass m 
on each side. 

By taking again advantage of the equal pressure condition, i.e., Tj^t) / X(t) = T B (t)/[L — 
X(t)] = 2T /L Eq.(|2} can be rewritten as 

' X-L/2 



jt + 4N kTorn / 1 | 1 ^ | 2Nm 



N V ttL V y/X y/L-XJ M 



X(L-X) 



X 2 = a(t). (3) 



Determining the correct expression for a(t) is a delicate task since the standard Langevin 
approach does not in general carry over to nonlinear dynamical systems jlfjj . as the one 
described by Eq.(jHJ). In order to be able to take advantage of the Langevin method, we 
linearize the above equation: a) by considering small displacements around the starting 
position of the piston, that is \(X — L/2)\/(L/2) << 1 and b) by approximating the square 
X 2 of the piston velocity by the average thermal value kT /M ; both hypotheses have to be 
proved consistent a posteriori. More precisely, by resorting to assumption b) and by using 
Langevin's method , it is possible to identify two characteristic times, t t h and t as » t t h 
(see Eq. (jlUJl ). The first one represents the temporal scale over which the piston velocity 
thermalizes, while the second one is the characteristic time over which the large piston- 
position fluctuations occur. The relation t th « t as justifies a posteriori the application of 
the fluctuation-dissipation theorem. 

After introducing the new variable x(t) = X(t) — L/2, the linearized version of Eq.fjHJ) 
reads 

vq = {2kTo/rn) l l 2 being the most probable velocity of the corresponding gas Maxwellian 



distribution 
theorem 



At this point, a straightforward application of the fluctuation-dissipation 



17( yields 

(a(t)a(t')) = ^v*(^ 2 S(t-n (5) 

According to the above results, the variable x(t), describing the evolution of the state 
of our system, obeys an equation formally similar to that of a "Brownian particle" in one 



dimension. However, the potential term and the viscous drag have naturally emerged out of 
the internal dynamics of the system and do not correspond to an external active force or to 
a phenomenological friction term as in the case of Brownian motors. If we define 

AN vq m 
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(6) 



L M ' L 2 \MJ 

Eq. (0J) is formally identical to that describing the Brownian motion of a harmonically-bound 
particle of mass M, that is 

x + fix + uj 2 x = a{t) (7) 

where a(t) is the Langevin acceleration, a case which has been thoroughly investigated in 
the literature 



III. THE PISTON AS A HARMONICALLY-BOUND BROWNIAN PARTICLE 



The solution of the stochastic differential equation (JJJ) 



governing the piston motion is 



181 ] for the Brownian motion of a 



readily obtained in terms of the general results provided in 
harmonically-bound particle. Equation (J7j) , together with Eq.(|5j). describe an "overdamped r 
case (in fact, (3/uj = ^Jl6/ny/N » 1, since N is typically a large number). For this situation, 
one has 



(x(t)) 
<«(*)> 



x e' 



u e 



-pt/2 
-pt/2 



( cosh M + A sinh + M e -W sinh M 
( C osh ^ - #- sinh - ^t e -Pt/2 ginh M 



(x(t) 2 ) = (x(t)) 2 + ^ [l - (2| sinh 2 f> + | sinh fat + 1 
(u(t) 2 ) = (u(t)) 2 + *g> [l - e"^ (2| sinh 2 f - | sinh^t + 1 

where tt(i) = dx/dt , xo = x(t = 0), wo — = 0) and 

fa = y//3 2 - Au 2 = (3- 2u 2 /(3. 



(8) 



(9) 



In our case xq = 0, uq = 0, so that, as expected, (x) = (u) = at all times (see Eqs. ©1,2)- 
By taking the temporal asymptotic limit of the above equations, it is possible to recognize 
the existence of two characteristic time scales tth and t as , such that tth « tas, given by 



tth — OR! 



2/3' 
2oj 2 • 



(10) 



7 



They respectively represent the time over which the piston velocity thermalizes, i.e., its 
mean-square attains the value 

(i\s = ^, (11) 

and the time over which the mean-square value of the piston position reaches the asymptotic 
value 

= ^ (12) 

By using Eq.© we have 



1 L M 1 L 

iV> (13) 



vq m 2y / 7rv 
and 

tth = ^t as . (14) 

Since iV >> 1, the asymptotic time t as is much larger than the thermalization time t t h, 
which justifies the replacement in Eq.fjSJ) of X 2 with its thermal value given by Eq. (jllj) . It 
is worth noting that the difference between the time scales of (x 2 ) and (x 2 ) (i.e. t as » t t h) 
is associated with the presence of the sign "plus" and "minus" in the r.h.s. of Eq. ©3 
and ©4, respectively. In fact, it is the minus sign in Eq. ©4 which is responsible for the 
cancellation of the slowly-decaying asymptotic terms in (u 2 ), while these terms, present in 
Eq. (jSJ) 3 , dictate the long time behavior of {x 2 ) as . 

The asymptotic mean-square displacement (x 2 ) as of the piston from its central equilibrium 
position x = is, according to Eqs. (jS)) and ©, 

2 \ V f L \ 2 1 M (L^ 2 



<^>- = 2 UJ = {-2) • < 15 > 

In the limit fi « 1 (that is, small piston mass with respect to gas mass), the root-mean 
square displacement (x 2 ) l J 2 is much smaller than L/2, so that both assumptions justifying 
the linearization of Eq.((7J), and thus the application of Langevin's approach, are verified. 
These fluctuations can be compared with those pertaining to a diathermal piston, associated 
with a system identical with the one described above but for the presence of a thermally 
conducting piston instead of an adiabatic one. To this end, let us assume the piston in Fig.l 
to be a good heat conductor, so that both sections possess the common constant temperature 
T . Any displacement x of the piston from the central position results in an unbalance p(x) 
between the pressures pa(x) and Pb{x) in the two sections, i.e., 

a / \ f , f , NkT NkT 

Ap(x) = p A {x) - p B {x) = A(L/2 + x) ~ A(L/2 _ x y (16) 



where A is the transverse section area, so that the force F(x) exerted by the gas on the 
piston reads 

Fix) - NkT ° - NkT ° - 2NkT \ (17) 
{ > L/2 + x L/2-x~ (L/2) 2 ' 1 ' 

having assumed \x\ << L/2, as verified a posteriori. Therefore, the piston feels both the 
thermal bath at temperature To and the harmonic potential 

= ^\ (18) 

and, as a consequence, its position probability-density P(x) reads 

N 



, * N 

p(x) = V^w exp 



:X 2 



(19) 



The associated mean square value (x ) is 



W'Mif' (20) 

smaller than the corresponding mean square displacement of the adiabatic piston (see 
Eq.(|15jl) by a factor m/M. Note that Eq. (|2()|) agrees with the standard result obtained 
in the frame of equilibrium-thermodynamic microscopic fluctuations. 

Our derivation clarifies the basic difference between adiabatic and diathermal situation. 
In the first case, the restoring force acting on the piston is given by —(8NkTom/L 2 M)x (see 
Eq.(jlJ)), while, in the second it reads —(8NkT /L 2 )x (see Eq. ()17|) ). so that the restoring force 
acting on the adiabatic piston is m/M times smaller than that acting on the diathermal 
one. We stress that the anomalously large fluctuations of the adiabatic piston are not 
limited to the macroscopic regime, but are also present in the mesoscopic regime, provided 
the number of molecules N is large enough, as in our case, to justify the use of ordinary 
concepts of pressure and temperature. This accounts for the peculiarity of the adiabatic 
piston system. The large random displacements of the piston are not trivially related to 
the system dimensions, but to its adiabatic nature: the diathermal piston exhibits a mean- 
square-value (x 2 ) = (m/M)(x 2 ) as (see Eq. (J2()jl and Eq. (jl5)l ) smaller by a factor m/M than 
the corresponding factor for the adiabatic piston. 



IV. VARIATION OF THE ENTROPY 



The peculiar behavior of the adiabatic-piston system discussed above has a natural coun- 
terpart in its entropy variations. In order to clarify this point, we note that the entropy 
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behavior of our system is not uniquely defined. First, let us consider an ensemble of identical 
systems, in each of which the piston is let indefinitely free to move starting from the central 
position x = at time t = 0. After a time interval reasonably larger than t as , the probability 
distribution of the piston position will stay unchanged. Accordingly, it is possible to define 
an entropy of the system (gas+piston) ensemble whose value turns out to be larger than the 
initial one pertaining to the one in which all pistons are fixed in the central position. The 
above entropy does not allow for a simple statistical interpretation like Boltzmann's one. 
The latter requires a single system in thermal equilibrium, and this is not the case for the 
wandering adiabatic piston. Thus, to give meaning to Boltzmann's entropy, we consider a 
single system and stop the piston either at a given time of the order of t as or at a given 
prefixed position of the order (x 2 ) l J 2 . The two procedures are conceptually quite different. 
In fact, halting the piston at a given time requires a device comparable in size with the 
piston-position fluctuations it is trying to harness: in other words, the stopping mechanism 
seems to behave like a demon device while it rectifies the large piston fluctuations, thus 
forbidding any violation of the Second Law. Similar arguments go back to Smoluchowsky 
(see ref. JjO]), and an interesting example is connected with the Feynman's pawl and ratchet 
device. |4| Vice versa, in principle, if the halting device is placed at a position x chosen a 
priori (e.g., x = +{x 2 )li 2 ), the piston stopping operation takes place without the interven- 
tion of relevant statistical fluctuations, i.e., without any demon mechanism. Of course, due 
to the random nature of the piston wandering, we do not know when the halting process will 
take place. However, if we wait long enough (e.g., for a time of the order of a few t as ), there 
is a high likelihood that the piston has come to a halt at the prefixed position. The above 
argument highlights a peculiar feature of the adiabatic piston system. Unlike other systems 
in which unusually large thermal fluctuations cannot be rectified without introducing com- 
parable external fluctuations, the piston wandering can be conveniently exploited. If the 
piston is stopped at a prefixed position, the Boltzmann gas entropy significantly decreases 
without any entropy increase of the environment. To clarify this point, we consider the 
entropy change AS" undergone by the system when passing from the equilibrium state in 
which the piston is held in the central position to the final one in which the piston is stopped 
at the position x = x. Recalling that the pressure does not change during the process (see 
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Sect. II), we have 

AO AO AO , L/2+x L/2-x 
Ab = Ada + Ad # = "ftCp In — ^—^ h nc D m — ^-^ — = nc„ In 



L/2 



L/2 



W2, 



(21) 



where c p is the molar heat at constant pressure and n the common number of molecules 
in sections A and B. By recalling Eq.(JT5J) and that x 2 /(L/2) 2 = (x 2 ) as /(L/2) 2 is a small 
number, Eq. (f2~T]) yields 



AS 



nc p 



In 



1 



x 



L/2 



nc„ 



x 



L/2, 



(22) 



or, by comparing it with the standard entropy fluctuations AS t h = (kncp) 1 ^ 2 |2jJ 

AS 



AS 



(23) 



This corresponds, whenever /iviV >> 1, i.e. M/m » y/N, to a sensible negative entropy 
decrease, and constitutes a violation of the Second Law, consistent with the perpetuum 
mobile behavior of the system. In particular, if the system is embedded in a thermal bath 
at temperature To, it is possible to devise a cyclic process by means of which the work 



W £ kT Nfi 



(24) 



can be extracted from the environment (see Appendix). Actually, according to the accumu- 
lated experimental evidence confirming the validity of the Second Law, we do not expect 
these results to hold true in standard macroscopic situations. In fact, our model imposes 
strict limitations on the spatial scale over which our results can be trusted. In order to assess 
these limitations, we refer to the specific case of a gas at standard conditions of temperature 
and pressure. In this case, expressing hereafter L in microns, and assuming the volume of 
the piston to be L 3 , N = 3 x 10 7 L 3 it follows from Eq.flHI that 

AS 



AS 



-5 x 10 3 L 3/ V 



th 



(25) 



while from Eq. (jl3|) we have, taking v o = 4 x 10 s microns/sec (molecular oxygen), 



t as = 5 x 10 2 /iL 4 sec. 



(26) 



We now observe that, if we reasonably require the work W extracted in a cycle to scale 
with the system dimensions, i.e., with N, the factor [i has to be kept constant when varying 
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FIG. 2: Normalized asymptotic standard deviation rj of the piston position as a function of u 
and N = 3 x 10 (averaged over 1000 realizations). The continuous line represents the predicted 
behavior in the linearized regime (see Eq. Q15JI ). 

the system dimensions. Therefore, Eq. (|2l)j) reveals an extremely sensitive fourth-power 
dependence of t as on the linear dimension of the system, which severely limits the appli- 
cability of our model to values of L up to a few microns. In fact, beyond this mesoscopic 
scale, t as becomes so large as to invalidate our results. As an example, for fi = 0.01 and 
L = lcm, one gets t as = 5 x 10 12 sec, that is about 1000 centuries! Conversely, by taking 
L = lfim, we get for t as a reasonable value of 5 x 10~ 4 sec, while AS/AS t h = —50 is still 
quite large. It is important to stress that this specific mesoscopic regime, around l[im, has 
emerged spontaneously from the adiabatic-piston problem, without the introduction of any 
ad hoc parameter. The Second Law, stating that the entropy of an isolated system cannot 
decrease, is actually violated in the specific example provided above. The associated spatial 
scale appears to set a borderline between the macroscopic realm and the mesoscopic one. 

V. BEYOND THE LINEAR REGIME 

The above approach has allowed us to deal with the situation /i < 1. In order to have an 
insight into the behavior of our process in the more general case \x ^$ 1, we assume Langevin's 
approach to be approximately valid also in this moderately nonlinear regime, and use the 
nonlinear Eq.© with the same stochastic acceleration a(t) worked out in the linear case. 
After introducing the dimensionless units £ = X/L and r = t/t Q , where t Q = A\/2t th , Eq.(jSJ) 
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FIG. 3: Time evolution of the r.m.s. of the normalized piston velocity u = d^/dr as a function of 
t' = 10 _4 t/t0) for = 0.5 and N = 3 x 10 4 , averaged over 1000 realizations. 



reads 



d 2 £ 



+ 



+ 



1 I 



1 



dr 



aa(r), 



(27) 



where a(r) is a unitary-power white noise process and a 2 = (7r/2 v / 2)( y u/A^). This last 
equation can be numerically integrated by adopting a second-order leap-frog algorithm as 



the one developed in reference 



22 



In particular, we can evaluate the asymptotic value of the mean-square root deviation 
r\ = ((£ — 1/2) 2 ) 1//2 as a function of [i and, for a given value of [i , the time evolution of 
{(d^/dr) 2 ). Fig.2 shows that the asymptotic mean-square root displacement of the piston 
increases as ^JJi for \i < 0.3 , in fairly good agreement with the linearized theory (see 
Eq.([15|)). while it tends to saturate for larger values of p. Fig. 3 confirms that the piston 
velocity attains its thermal value in a time much shorter than the asymptotic time t as , 
a necessary condition for the application of Langevin's approach. As far as the entropy 
change AS is concerned, the linear dependence on x 2 = (x 2 ) as predicted by Eq. (j22)l implies 
as well its saturation for increasing /x . The above results show that, for fixed N, both piston 
asymptotic displacement and entropy decrease saturate with increasing piston mass M. 
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VI. CONCLUSIONS 



In the present paper, we have examined the dynamic evolution of the adiabatic piston 
system starting from both mechanical and thermal equilibrium. Obviously, in this case 
(x) = and (x) = for symmetry reasons, so that the significant quantity is {x 2 (t)) l l 2 . 
This situation can be analyzed by generalizing the deterministic kinetic model developed in 
js] with the introduction of Langevin's force. The resulting equation for x(t) turns out to 
be, whenever the mass M of the piston is small compared to the gas mass M g , completely 
identical to that describing the Brownian motion of an harmonically-bound particle of mass 
M. This allows us to provide an analytic expression for {x 2 (t)) l l 2 , which can be a size- 
able fraction of the cylinder length L (see Eq. (|15jl ). thus implying that the piston never 
reaches an equilibrium position but keeps oscillating around its initial position x — 0. This 
peculiar behavior is restricted to systems possessing specific spatial dimensions pertaining 
to a mesoscopic realm, where the laws of thermodynamics are not necessarily valid. The 
fact that the piston never stops is strictly related to the finite length of the cylinder: if we 
allow the length L of the cylinder to go to infinity (as assumed, for example, in jjj), the 
ratio M/M g goes to zero and so does {x(t) 2 ) l l 2 /(L/2) (see Eq.(jl5jl). At the same time, t as 
becomes extremely large, so that, in this limit, the piston practically does not move. In view 
of this, the violation of the Second Law does not appear that startling since it concerns a 
regime outside the macroscopic thermodynamic limit. However, it can be interpreted as the 
signature of something remarkable happening at the borderline between the macroscopic 
and mesoscopic realms, signaling that a new way of looking at standard thermodynamic 
concepts may be needed if we wish to continue to apply them outside the boundaries within 
which they were first introduced. 

In order to have an insight into the regime of mild nonlinearity occurring when fi = 
M/Mg approaches unity, we have numerically solved the stochastic Eq. (j27|) with the standard 
initial conditions adopted in our model. The results of the linearized approach are in good 
agreement with those associated with the numerical solution of Eq. (|27j) . 

The same problem of time evolution and approach to equilibrium of the system has been 
also considered by using molecular dynamics simulations. [lcl llll . I12I lHI ] These typically 
involve a considerable number of point particles (around 500) which model the gas inside 
the cylinder and are separated by a frictionless movable piston, without internal degrees 
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FIG. 4: Cyclic process through which work can be extracted from a thermal bath at temperature 
T . 

of freedom, against which they undergo perfectly elastic collisions. Some general features 
emerge from the numerical analysis, which reveal, consistently with our results, a very slow 
approach to the final equilibrium state. 



In particular, Ref |12j presents the molecular dynamics simulations of a system evolv- 
ing from an initial state in which the piston is fixed and the number of particles and the 
temperatures in both sections are equal, which is precisely the case we are dealing with 
in the present paper. Our analytic approach predicts a relaxation time (see Eq. (|lHjl ) 
which, in the normalized units adopted in (vq — 2, m — 1 and L = 60), reads 



t as = (30/v/7r)(M/m) = 17 M. This is in remarkably 



good agreement with the molecular dy- 



namics simulations reported in [12] (see Fig. 5 of Ref. 12 1) where, for M < 100 (corresponding 
to n = M/N = 100/250 = 0.4), t as = 16M. 

In conclusion, our simple kinetic analytic description has revealed that Callen's adiabatic 
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piston can exhibit, when starting from a particular mechanical and thermal equilibrium 
state, large fluctuations, a very intriguing feature which presents the characteristics of per- 
petuum mobile. This occurs for very specific spatial dimensions of the system (around 1/im), 
pertaining to the mesoscopic regime. Whether this is actually challenging the limits of va- 
lidity of the Second Law seems a legitimate question, whose answer may shed some light on 
the understanding of small-scale nonequilibrium devices. 

We wish to thank Noel Corngold for his constant encouragement and many useful com- 
ments. 
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Appendix 

If our system is embedded in a thermal bath at temperature To and pressure Pq, the 
following process, sketched in Fig. 0], can be used for extracting work from the environment. 
In the first panel, the system is in the initial state. The second panel shows the system after 
the entropy decrease AS has occurred. At this point, we separate the two sections (panel 
3) and drive the two gases back to the initial temperature To by means of two reversible 
adiabatic processes, the relation between the new volumes V{ and V 2 ' and the previous 
volumes V\ ,V 2 and temperatures Ti,T 2 being 

V{ = V 1 (T 1 /T )^, V> = V 2 (T 2 /T )^t. (28) 

This is followed (panel 4) by two reversible isothermal processes in which the two sections 
are in contact with the thermal bath at temperature T , which take the two sections back 
to the initial temperature T and volume Vo/2 (frame 5). The total work extracted in the 
process is 

W = W 1 + W 2 = nRT In (-£ + nRT ln(-£) = nRT In ol 0> . (29) 

V i V 2 V^iT.T,) — 

Since V\ + V 2 = 2V and T x + T 2 = 2T , we obviously have Tq > T 1 T 2 and Vq 2 > ViV 2 , so that 
a work W larger than zero has been obtained at the expenses of the heat Q = W extracted 
by a single source (the thermal bath). It can be easily checked that W = —To AS* = kToN/i, 
that is each gas molecule contributes with a fraction /xfcTo/2 to the work extracted in the 
process. 
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